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O ; ABSTRACT 

> ■ We analyze the central dark-matter (DM) content of ~ 4, 500 massive (M* >10 10 M@), 

\ low-redshift (z < 0.1), early-type galaxies (ETGs), with high-quality ugrizYJHK 

photometry and optical spectroscopy from SDSS and UKIDSS. We estimate the "cen- 
tral" fraction of DM within the if-band effective radius, R c s, using spherically sym- 
metric isotropic galaxy models. We discuss the role of systematics in stellar mass 
estimates, dynamical modelling, and velocity dispersion anisotropy. The main results 
of the present work are the following: (1) DM fractions increase systematically with 
. both structural parameters (i.e. i? c ff , and Sersic index, n) and mass proxies (central 

■ velocity dispersion, stellar and dynamical mass), as in previous studies, and decrease 

with central stellar density. 2) All correlations involving DM fractions are caused 
by two fundamental ones with galaxy effective radius and central velocity dispersion. 
These correlations are independent of each other, so that ETGs populate a central-DM 
plane (DMP), i.e. a correlation among fraction of total-to-stcllar mass, effective radius, 
and velocity dispersion, whose scatter along the total-to-stellar mass axis amounts to 
~ 0.15 dex. (3) In general, under the assumption of an isothermal or a constant M/L 
profile for the total mass distribution, a Chabrier IMF is favoured with respect to 
r j2j ■ a bottom-heavier Salpcter IMF, as the latter produces negative (i.e. unphysical) DM 

& \ fractions for more than 50% of the galaxies in our sample. For a Chabrier IMF, the DM 

estimates agree with ACDM toy-galaxy models based on contracted DM-halo density 
profiles. We also find agreement with predictions from hydrodynamical simulations. 
(4) The central DM content of ETGs does not depend significantly on the environment 
where galaxies reside, with group and field ETGs having similar DM trends. 

Key words: dark matter - galaxies : evolution - galaxies : galaxies : general - 
galaxies : elliptical and lenticular, cD. 
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1 INTRODUCTION 



the 



In 

(Slo an 



last decade, 
Digital Sky 



large ar ea surveys like 
S urvey; [Abazaiian et al.l 



SPSS 
20031 ; 



lAdelman-McCarthv et~aH 120081 ; lAbazaiian et al.1 l2009t ). 

have provided high quality data contributing significantly 
to our understanding of galaxy properties and scaling 
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relations. Nevertheless, galaxy formation remains one of the 
outstanding questions of modern astrophysics. As the most 
massive stellar systems in the nearby universe, early-type 
galaxies (ETGs) have a special role in providing the under- 
pinnings for a consistent galaxy formation picture. They 
form a relatively homogeneous class of objects, dominated 
by an old stellar population, a small fraction of cold gas 
and low- levels of ongoing star formation (SF). These char- 
acteristics make them a potentially powerful tool to trace 
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the evolution of cosmic structures back through the cosmic 
epochs. The uniformity in the ETG properties involves tight 
correlations of quantities like effective radius, Rett, surface 
brightness measured within this radius, and central velocity 
dispersion, op, which merge into the so-called Fundamenta l 
Plane (FP; IDiorgovski fc Dav"islll987l ; IDressler et al.lll987l ). 
An observed deviation between the FP coefficients and 
those expected from the virial theorem has been interpreted, 
among other possibilities, as a variation of the tot al M/L 
with galaxy luminosity /mass (|Dressler et al.1 1 19871 ') which 



Ciotti et al.1 1 19961; iBusarello et al.1 1 19971; iGraham & Colless 


1997 


Prugniel & Simienl 


19971: Truiillo. Burkert & Bell 


2004 


Cappellari et al. 2006 


; D'Onofrio et all 2006 


; Graves 


2009 


Tortora et al.1 12009. T+09; lLa Barbera et al 


l2010bh. 



Therefore, measuring the DM content in an independent 
way is of paramount importance for studying the ETG 
scaling relations as well as the overall process of galaxy 
formation and evolution. Recently, the DM content in the 
central regions of ETGs (typically wi thin lR B s) has been 
analy z ed using both local sa m ples ([padmanabhan et al 
12004 : ICappellari et al.1 120061; Hyde fc Bernardil 



T+09; 
after 



Napolitano. Romanowskv fc Tortora I 2010l . here- 
NRT10) and intermediate-r e dshift gravi t ational 



lens galaxies ([Cardone et al. 20091; lAuger et al.1 l2010bl: 



ICardone fc Tortora |201 



. iTortora et al.1 l20ic[ hereafter 
T+10; IFaure et al.1 1201 ll ; IMore et al.1 l201lf ). Independent 
of the model used to describe the mass distribution in 
a galaxy, several studies have found that DM fractions 
within J? e ff increase with galaxy luminosity, stellar mass, 



size, a nd velocity dispe r sion ( e.g., Ferreras. Saha. Williams! 



Forbes et al.1 120081 ; T+09; 



120051 ; ICappellari et"ai1 12006.. , 
NRT10; lAuger et al.ll2010bl ; iLeier et al.ll201ll; T+10), while 
a diffe rent conclusion was drawn by Truiillo. Burkert fc Belli 
(2004), base d on a constant- M/L mass model (see also 
iGrillol l2010h . The correlation of DM density with stel- 
lar mass and R e ft points to DM profiles being cuspy 
(|Thomas et al.1 |2009| ; T+09; NRT10; T+10). However, 
even in this case, inconsistent findings have been reported. 
Using rotation curves of spiral galaxies, and mass models 
of individual dwarf and spiral galaxies, as well as the 
weak lensi ng signal of ell i ptical s and spirals, iDonato et al.l 
(120091 ) and lGentile et al. I (|2009l ) found that the central DM 
column density is constant over twelve orders of magnitude 
in luminosity. On the contrary, NRT10 showed that, on 
average, the projected central density of nearby ETGs is 
systematically higher than that of spiral and dwarf galaxies, 
implying an increase of DM dens ity with halo ma s s. Th e 
same conclusion was reached by iBovarskv et al. I (|2009l ). 
using data for different galaxy types and groups/clusters of 
galaxies. Part of the above controversies may be associated 
with sample selection issues, with differences in the way 
galaxy parameters are measured, as well as with different 
assumptions a bout the halo models, I MF, and adopted 
fiducial radius (jCardone fc Tortorall201ol ). In this paper, we 
study how the central DM content of ETGs correlates with 
galaxy properties, for a large and homogeneous sample 
with a wealth of photometric and spectroscopic data 
available (|La Barbera et al"1l2010al . Paper I of the series). In 
previous papers, we analyzed the FP rela tion of ETGs and 
its de pendence on galaxy environment (lLa Barbera et al.l 
l2010bl . Paper II; lLa Barbera et al.1 l2010d . Paper HI); 



and the correlation of in ternal colour gradients o f ETGs 
with galaxy properties (|La Barbera et ah! l2010bl . Paper 
IV). Here, we investigate the DM content of ETGs in 
terms of both structural parameters and various mass 
proxies, contrasting the average trends with predictions of 
toy-models and cosmological simulations. We show, for the 
first time, how the correlations between the DM content 
and galaxy properties depend on the environment where 
these systems reside. Our study resembles previous ones 
in its use of SDSS data to stu dy the connections b e tween 

stellar and dynami cal ma ss ([Padmanabhan et al.1 2004 • 

Shankar fc Bernardil 120091 ; [G raves fc Faberl \2oW IGrillol 
20101) . but goes beyond these by incorporating additional 
photometric data as well as environmental information, 
and may be considered the definitive study of the central 
DM content of a large sample of bright ETGs in the SDSS 
era. The outline of the paper is as follows. In Section [2] we 
describe the sample and the data analysis as well as stellar 
and dynamical mass calculations. In Section [3] central DM 
fraction and density are analyzed as a function of galaxy 
mass and structural parameters. In Section U we discuss 
the main drivers of the DM content and the "DM plane" 
of ETGs, and in Section [5] the implications for the FP are 
discussed. Section |(J] deals with galaxy environment, while 
Section [7] compares the correlations involving DM content 
with expectations from cosmological and toy-galaxy models. 
Conclusions are presented in Section [5] Throughout the 
paper, we adopt Hq = 75kms _1 Mpc - . 



2 DATA SAMPLE AND ANALYSIS 
2.1 Sample 

A volume-limited sample of 39, 993 "bright" ( 01 A/ r < -20) 
ETGs, in the redshift range of 0.05 to 0.095, with available 
ugriz photometry and optical spectroscopy from SDSS-DR6 
is used in this work and in previous papers of the SPI- 
DER (Spheroid's Panchromatic Investigation in Different 
Environmental Regimes) project. 5, 080 galaxies also have 
YJHK photometry from DR2 of UKIDSS-LAS (see Pa- 
per I). ETGs are defined as bulge dominated systems, with 
passive spectra in the ir centres (i.e. within t he SDSS fibre 
apertures). Following iBernardi et all (|2003al ). from an op- 
erational viewpoint, ETGs are those systems with eClass 
< and f racDevr > 0.8, where the SDSS spectroscopic pa- 
rameter eClass gives the spectral type of a galaxy, while 
the SDSS photometric parameter fracDevr measures the 
fraction of galaxy light that is better fitted by a de Vau- 
couleurs (rather than an exponential) lawQ. All galaxies 
have central velocity dispersion, a, from SDSS-DR6, in the 
range of 70 to 420 km s - . In all wavebands, galaxy struc- 
tural parameters - i.e. the effective radius, R e g, the mean 
surface brightness within that radius, < (x > e , and the 



1 Notice that the eClass and fracDevr selections are very effec- 
tive to remove late- type systems (see Paper I), but do not allow 
a clear separation of E and SO galaxy types. Since SOs are flat- 
ter than ellipticals, we have carried out a test of restricting the 
analysis to objects with it-band axis-ratio q > 0.8. With this 
alternative cut that minimizes the fraction of SOs, we find no sig- 
nificant variation in the results presented throughout the paper. 
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Sersic inde x, n - have been homoge neously measured using 
2DPHOT \L& Barbera et al.ll2008ah , fitting galaxy images 
with seeing convolved two-dimensional Sersic models. The 
SPIDER sample is 95% complete at 1 M r = -20.32 or, at 
a stellar mass h'L = 3 x 1O 1O M for a lChabrier I (120011 ) IMF. 

For the present study, we select SPIDER ETGs with 
high quality structural parameters in the optical and Near- 
InfraRed (NIR) wavebands, according to the following cri- 
teria: the Sersic fit has x 2 < 2 in all wavebands; uncertainty 
on log R e g < 0.5 dex from g through K (see Paper IV for 
details on these thre sholds); available stellar mass estimates 
((Swindle et al.|[201ll . Paper V). The resulting sample con- 
sists of 4, 259 ETGs. As expected, all these galaxies reside 
on the red-sequence, with more than 99% having g — r J>0.5 
(within an aperture of 1 i? e ff), and a median g — r = 0.88. 
For this sample, the R e g spans the range ~ (0.5 — 40) kpc 
( ~ (1 — 15) kpc for more than 90% of the galaxies), with 
a median of 3.5 kpc. The median ratio between the SDSS 
fibre aperture, where velocity dispersions are measured, and 
the K-b&nd effective radius, R^/Reft, amounts to ~ 0.6, 
implying only a little "extrapolation" in our DM fraction 
estimates (see Sec. 12. 3[) . As shown in Paper I, significant 
differences are found when comparing 2DPHOT and SDSS 
structural parameters. Such differences arise from the use 
of Sersic (2DPHOT) rather than de Vaucouleurs (SDSS) 
models to fit the light and total mass distribution of ETGs; 
from different software to measure the structural parame- 
ters; and fro m the sky estimation bias tha t affects SDSS 
photometry (|Adelman- McCarthy et al]|2008l ). In Appendix 
A, we further illustrate this point by comparing correlations 
among structural parameters from different sources. The en- 
vironment of ETGs in the SPIDER sample is characterized 
by a friends-o f -frien ds catalog of 8, 083 groups, created as in 
iBerlind et afll|2006h . Here, we used Berlind's algorithm over 
a larger area, as the new SDSS-DR7 is now available (rather 
than DR3; see Paper III for details). A shifting gapper tech - 
nique is applied to this catalogue (see lLopes et al.ll2009al ). 
allowing galaxies to be classified as either group members 
(~ 46%), non-group members (hereafter "field" galaxies; 
~ 33%), or unclassified (~ 21%). We also separate group 
members into central and satellite galaxies, where the cen- 
tral galaxy o f a given group i s the object with the largest 
stellar mass jYanget al.ll2007l ). 

2.2 Stellar mass estimates 

Stellar masses are derive d by fitting synthetic stella r pop- 
ulation (SP) models from iBruzual fc Charlotl (|2003L BC03 
hereafter) to th e optical+NIR pho tometry, using the soft- 
ware LePhare jllbert et al.l 120061 ). The observed galaxy 
fluxes determine the normalization of the best-fit template, 
which then gives the stellar mass, M* (initially computed 
within the SDSS Kron aperture, which we correct using 
Sersic models to total mass). We adopt a set of SP models 
with different star-formation e- folding times (r ^ 15 Gyr), 
internal reddenings (E(B — V) ^ 0.5), and metallic ities 
(0.2 ^ Z/Zq < 3), assuming a ICardelli et alj (|l989l ) ex- 
tinction law and Chabrier IMF. 

We refer the reader to Paper V for all details about 
the fitting procedure and the accuracy of stellar mass es- 
timates using different theoretical and empirical assump- 
tions. We provide here only a summary of statistical and 



systematic uncertainties on stellar mass estimates. In paper 
V, we compared stellar mass estimates obtained with dif- 
ferent st ellar population models, e.g. BC 03 and PEGASE.2 
models |Fioc fc Rocca-Volmerangelll997l ). different internal 
extinction laws, and different combinations of wavebands 
adopted in the SED fitting procedure. Also, we performed 
comparisons with spectroscopically derived stellar masses 
(the two sets of estimates being in excellent agreement, with 
a median difference of 0.03 dex). 

In general, the comparison of different sets of stellar 
mass estimates shows that the scatter in M* (i.e. its un- 
certainty) ranges from 0.05 to 0.15 dex. As in other studies 
(see Paper V and references therein) systematic uncertain- 
ties on stellar mass were found to play a major role, with 
variations in the IMF and extinction law yielding systematic 
biases on the mass of nearly a factor of 2. Notice though, 
that despite the age-metallicity degeneracies in photometric 
data, these conspire to keep the stellar mass-to-light ratio 
(and hence the stellar masses) relatively well constrained. A 
de-projected Sersic law in the AT-band is used to describe 
the density profile of the stellar component. 

2.3 Dynamical mass estimates 

The dynamical mass within a given galaxy radius, r, is usu- 
ally computed as: 



Mdyn = 



Ka 2 r 



(1) 



where G is the gravitational constant, a is the galaxy ve- 
locity dispersion, and K is a pressure correction term (e.g. 
Padmanabhan et al.ll2004l ; lEke et alJkoO^ ICappellari et al.1 



20061 : T+09). The K depends on several factors, like the ra- 



dius wherein Mdyn is computed, the aperture used to mea- 
sure <t, the viewing angle of the system, its orbital structure, 
luminosity profile, and how DM is distributed. Here, we fol- 
low the approach used in T+09, where instead of adopting 
some approximation for K in Eq. {TJ, we model each indi- 
vidual galaxy directly using the J eans equations to estimate 
Mdyn within r = 1 7? e ff (see also ICardone fc Tortoral |2010| ; 
T+10). The models also require assumptions about the mass 
profile in order to extrapolate to R e s from the more central 
(on average) a measurements. Here our approach is to adopt 
two types of mass profile that bracket a range of possibili- 
ties. Each of these describes the total mass profile with one 
free parameter: 

• SIS. The Singular Isothermal Sphere, where M(r) oc 
o"s Is r, is the prototype of a galaxy profile producing a flat 
rotation curve. Despite its simplicity, this reproduces quite 
well the total mass profile in the mos t massive ETGs (e.; 
lKochanek1ll99ll : lKoopmans et af 1 l2006l : lGavazzi et al.ll2007 
and in particular the massive halo present in their outer re- 
gions , consistently with virial mas s estimates ([Benson et al.l 
200d: iMarinoni fc Hudson! 120021 : iNapolitano et al.1 1200a 



van den Bosch et al.ll2007T ). Since a galaxy halo is expected 



to be truncated when it enters a group/cluster of galaxies, 
we analyze the impact of the environment on the central 
DM fractions using also a truncated SIS model (tr-SIS), 
where the density profile is truncated at a radius rt- For 
each of the group ETGs, the value of rt is computed from the 
projected distance between the ETG and the centre of th e 
parent group, following the recipes of iGhigna et all (| 19981 ). 
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Figure 1. Panel (a). Distributions of K for SIS and const-M/L mass models (see Eq. l[TJ). Panels (b-e). Median value of K as a function 
of <7 (b), Mi, (c), R e ff (d), and n (e). Median values, with error bars showing 2 5—75 per cen t scatt er. Red and blue colours refer to SIS 
and const-M/L models, respectively. The black curve in panel (e) is taken from lBertin et all (120021 s ) . 



• Constant M/L. Some studies of ordinary ETGs suggest 
that the mass profile falls off roughly as steeply as the stars, 
at small to intermediate radii, which may imply either little 
DM in their centres or DM profiles that mimic the stars 
dRomanowskv et al 1l2003l : iDouelas et al1l2007l ; iDeason et al.l 
|2012| ). The mass profile is given by M(r) = To L(r), where 
To is the typical mass-to-light ratio of the sys tem and L(r ) 
is the luminosity profile, modeled here by the ISersic] (|l96&f ) 
law. Hence, the mass distribution is assumed to follow the 
profile of the stellar component, whose shape (i.e. the Sersic 
n), changes from galaxy to galaxy, reflecting the structural 
non-homology of the light profiles of ETGs. 

Operationally, for each galaxy, given a mass profile 
model, we compute the projected mass-weighted velocity 
dispersion of the model, <r mode i, within the SDSS fibre aper- 
ture. The a mo dei is then matched to the observed central 
velocity dispersion of the galaxy, providing the correspond- 
ing best-fit parameter, crsis (To) for the SIS (const-M/L) 
model, and the mass profile Md yll (r) (see Appendix [Bl for 
details) . 



3 DYNAMICAL MASSES AND CENTRAL 
DARK MATTER 

We characterize the DM content of an ETG by computing (i) 
the de-projected DM fraction, /dm(^) = 1 — M-,,{r) / M Ayn {r) 
or the total-to-stellar mass ratio Md y n(r)/M*(r),where M* 
and Mdyn are the stellar and dynamical mass □ within a 
given de-projected radius, r; and (ii) the de-projected aver- 
age DM density, (p DM ) = (M dyn (r) - M*(r))/(f 7rr 3 ), giv- 
ing more direct information on the DM content of a galaxy. 



The M*(r) is estimated by de-projecting the if-band light 
profile of each galaxy (assuming spherical symmetry), and 
normalizing this de-projected profile to the total M* esti- 
mate, from Paper V. Hereafter, we refer to /dm(>~ = Re) 
{{pDM{r = Re)}), i.e. the DM fraction (density) computed 
within a de-projected radius equal to the projected R e , as 
the "central" DM fraction (density). 



3.1 Dynamical masses and virial coefficients 

To provide a more straightforward way of estimating Md yn 
from velocity dispersion and effective radius (as done in pre- 
vious work), and of comparing mass estimates from different 
models, we also recast the computation of dynamical masses 
in terms of Eq. (TTJ . For each galaxy, we cor rect the observed 
cr to a n apertur e of 1 R e s , cr(-R e ff), followinglCappellari et al.l 
(|2006h (see also |j0rgensen et al.|[l995l . ll996l '). Then, we insert 
RcH, o~(Res), and the dynamical mass estimates obtained 
through the Jeans equations (see above) into Eq. (UJ) to ob- 
tain the corresponding virial coefficients, -Ksis and -K M / L , 
for the SIS and const-M/L models, respectively. The distri- 
butions of Ksis and -K M / L are displayed in the left-panel 
of Fig. [TJ while the right panels show the median values of 
-Ksis and .Km/l as a function of <r, M*, i? c ff, and n. The 
distributions of Ksis and K M / L have median values of 
Ksis = 2.571^7 and K m/l = 1.85±g;|| (where error bars 
are the 25-75 per cent scatter). The Km/l, is always smaller 
than -Ksis, reaching equality only for galaxies with a K- 
band Sersic n ~ 3.5 (see panel (e) of Fig. [TJ. This is con- 
sistent with the findings in T+09, where galaxy structural 
parameters were derived using a de Vaucouleurs rather than 



2 We implicitly assume that the dust and gas components give a 3 Using an aperture of i? fj/8, rather than li? e ft i gives i^siS 

negligible contribution to M^ yn . 2.14^^'^ and ifjyi/L = 1-54^q'|j. 
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a Sersic model. From Fig. Q] we see that Ksis is indepen- 
dent of a, AA, R D a, and n, while A M / L decreases with M*, 
Reg, and n, and is independent of a. The K — n trend is in 
fairly good agreement w i th th e best-fitting K — n relation 
obtained bv lBertin et al.l |2002l ) (see solid black line in panel 
(e)). Notice that using a constant Ksis rather than fitting 
each individual galaxy with the SIS model would introduce 
an uncertainty in dynamical mass equal to the scatter of the 
Ksis distribution in Fig. [T] ~ 4 — 7% (taking the lower and 
upper 25-75 per cent range, respectively). Similarly, using 
the K — n median trend to obtain A' M / L would imply an 
uncertainty of ~ 4-5% on Md yn - Modelling each individual 
galaxy has the advantage of avoiding this source of scatter, 
which could be important for discerning subtle differences 
in DM content among different subsamples (e.g. different 
environments, see Sec. [6|. 

For group ETGs, we have also computed K values by 
using the truncated SIS model, tr-SIS (see Sec. l2.3|) . finding 
them to be fully consistent with those obtained from the SIS 
model. In particular, only < 2% of group galaxies exhibit 
a difference of A log A' > 0.01 (~ 2.3%). This is due to 
the fact that the truncation radius, rt (see Sec. 12. 3[) . has a 
distribution peaked at ~ 50 kpc, with R c g being ~ 3.5 kpc, 
i.e. < r t . 

Throughout the present paper, the DM content of ETGs 
is estimated under the assumption that their stellar orbits 
are isotropic, which is incorrect at some level. Although a 
detailed analysis of anisotropic orbits is far from being triv- 
ial, and is certainly beyond the scope of the present work, 
we have estimated how anisotropy can affect our dynamical 
mass estimates. 

Detailed dynamical modeling efforts of ETGs have fo- 
cused extensively on their central regions (i.e. those we in- 
vestigate in this work), and found anisotropies to be fairly 
mild in general, typically in the range —0 -2 ^ ft ^ +0.3 
I Gerhard et al. hoOlMCappellari et alj|2007f ), where f3 = 1 — 
uJ/of quantifies the relative internal dispersions in the tan- 
gential and radial directions. Mild central aniso tropy is also 
predi cted from simulations of merger remnants |Dekel et al.l 
l2005h . 

Recent work has focused attention on the ability to 
constrain Mdyn(-Reff) in stellar systems, independently of 
anisotropy, if th e velocity dispersion can be measur ed over 
the whole galaxy l |Walker et alj|2009l ; lwoff et alj2010h . If the 
dispersion is measured only within an aperture of R c fs, but 
the dispersion profile is relatively constant with radius, it 
is possible that the mass can still be fairly well constrained 
(G. Mamon, priv. comm.). To check this, we have carried 
out test models as in Appendix B, but adopting constant 
anisotropy profiles, f3(r) = /?. 

For /3 = 0.1 (= 0.2), we find dynamical masses that 
are ~ 2% (~ 4%) smaller than in the isotropic case. In the 
case of extremely (and indeed unrealistic) radial (tangential) 
orbits, i.e. f3 — 1 fj3 = — 1), the masses are underestimated 
(overestimated) by ~ 30% (~ 10%). The effect is larger 
for the most massive and largest galaxies in our sample, 
for which we find that /3 = 0.1 (= 0.2) would make the 
masses smaller by ~ 3% (~ 6%), while in the case of /3 — 1 
(/3 = —1), the masses would be smaller (larger) by ~ 50% 
(~ 20%). Thus, for typically observed values of {S (~ 0.2), 
we do not expect any significant variation in the DM content 
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Figure 2. Dynamical mass within 1 R c ff as a function of stel- 
lar mass within the same radius. The black solid (dashed) line is 
the median SIS (const-M/L) trend, with the shaded region cor- 
responding to the 25-75 per cent scatter for the SIS case. The 
red lines are the median trends obtained by using only galaxies 
with /dm 0. A Chabrier IMF is assumed in the computation of 
AA (see Paper V). The solid blue line is the one-to-one relation 
-^dyn(^cff) = AA(Roff)' The arrows show the variation of M* 
when adopting a Salpetcr, rather than a Chabrier IMF. 

of ETGs as inferred under the assumption of isotropy, the 
effect being <^ 5%. 

3.2 Dark matter fractions 

Fig. [2] plots median dynamical masses as a function of stel- 
lar mass, M+(R e ff), for both SIS and const-M/L models. 
We find M dyn (R cS ) oc AA(A off ) a , with a > 1, i.e. the cen- 
tral DM content increases with galaxy mass. In particu- 
lar, we find «sis ~ 1-2 and om/l ~ 1-07 for the SIS and 
const-M/L models, respectively. Notice that on average, at 
fixed AA(A e fT), the A/ dyn (7? cf r) is larger than M*(i? ef f), and 
the fraction of ETGs with AA>Af dyn (i.e. unphysical DM 
fractions) is ~ 12% (24%) for the SIS (const-M/L). Using 
a Salpeter (rather than Chabrier) IMF would shift M* to 
larger values by ~ 0.25 dex (see the blue arrows in the Fig- 
ure) , making the median Mdyn only slightly larger than AA , 
and increasing the fraction of objects with AA>Afdy n up to 
~ 55% (78%) for the SIS (const-M/L) model. 

Taken at face value, this result implies that the bottom- 
heavy Salpeter IMF is disfavoured with respect to a Chabrier 
IMF, in agr eement with previous findings from stellar dy- 
namics fe.g..lCappellari et al.ll2006l; NRTIO). and fro m lens- 
ing l|Ferreras. Saha. Buries 20081; iFerreras et al]|2010h . How- 
ever, we warn the reader that (i) most of the objects with 
AA>A/d yn are from a region of parameter space (low AA) 
where our sample is more affected by incompleteness (see 
below); and (ii) our analysis relies on the assumption of a 
given model for the dynamical mass distribution in ETGs 
(i.e., either SIS or const-M/L). Indeed, if we relax the as- 
sumption of a universal IMF, our data would be consistent 
with eit her a Salpeter or Chabrier IMF a t high mass (T+09, 
NRT10, lAuger et all l2010al : iTreu et all l201Ch and with a 
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Chabrier IMF preferred at low mass ([Barnabe et al. 20 111; 
Sonnenfeld et all 2011 1_ Brewer et alj|2013 ; Cappellari et al.l 
20121 : button et al.ll2012l : see also Section EHl below) . 

Fig. [3] plots central DM fractions as a function of dif- 
ferent galaxy parameters, i.e. effective radius, Sersic index, 
velocity dispersion, stellar and dynamical mass, and cen- 
tral average de-projected stellar density, (p*), defined as 
(/£>*) = M ir (R e g)/(^nR e g 3 ). Black solid (dashed) lines show 
median trends for the SIS (const-M/L) model, while grey 
regions are the 25-75 per cent quantiles (i.e. the scatter) 
for the SIS model. Notice that, in general, the scatter goes 
from ~ 0.2 at low /dm (~ 0.2) up to ~ 0.1 (or smaller) at 
high /dm (~ 0.8). This is likely caused by the statistical un- 
certainties on galaxy parameters, such as, in particular, the 
Mi,. In fact, a typical statistical uncertainty of ~ 0.1 dex on 

(see Sec. l2.2[) propagates into an error of ~ 0.18 (~ 0.05) 
on /dm, for /dm~ 0.2 (/dm~ 0.8). 

For each galaxy parameter, we derive the completeness 
limit of the sample with respect to that parameter by looking 
at its correlation with the SDSS r-band Petrosian magnitude 
(which is the main parameter used to select SDSS spectro- 
scopic targets; see Papers II and IV for further details on 
this approach) . These completeness limits Q are marked by 
vertical dotted lines in Fig. Notice that because of the 
scatter around the correlations between /dm and the dif- 
ferent galaxy parameters, the range of /dm, as seen by the 
black curves in Fig. is different for different quantities, 
being larger for _R e ff , (p*) and M dyn (0 < /dm < 0.8), than 
for M^, a and n (0.2 <^ /dm <i 0.6). All median trends in 
the Figure are modeled by linear least-squares fits, adopt- 
ing log /dm as the dependent variable, and using only the 
portions of each diagram where the sample is complete 0. 
In these ranges, the log /dm trends are reasonably well de- 
scribed by the linear fits. The slopes of the best-fitted re- 
lations are reported in Table [1] where it should be kept in 
mind that the range of completeness corresponds to those 
galaxies which required the most extrapolation from the fi- 
bre aperture to -R e ff , so the /dm results are the most model 
dependent. We also tabulate the slopes obtained by using 
log A/dyn/A/*, rather than log /dm, as the dependent vari- 
able in the fits. This is relevant to connect DM fractions to 
the origin of scaling relations of ETGs (Sec.JSJ. Notice also, 
that due to scatter in DM content at a given point of the 



Table 1. Slopes of the correlation between /dm, M dyn /Mi, 
and (pdm) vs R c g, n, a, Mi,, (pi,) and M dyn , for the SIS and 
const-M/L models. We also show in parenthesis the slopes when 
A^dyn > Mi, . The fits are performed taking into account the com- 
pleteness limits. 







SIS 


const-M/L 


/DM - Rett 




0.26 ± 0.02 


0.27 ±0.05 






(0.24 ± 0.02) 


(0.18 ±0.04) 


/dm - n 




0.39 ± 0.11 


-0.60 ± 0.11 






(0.33 ± 0.07) 


(-0.22 ± 0.06) 


/dm — a 




0.71 ± 0.11 


1.17 ±0.18 






(0.62 ± 0.10) 


(0.84 ±0.15) 


/dm - Mi, 




0.24 ±0.04 


0.29 ±0.06 






(0.15 ± 0.03) 


(0.07 ±0.05) 


/DM — (p*) 




_n 11 -4- D 01 


_n 11 m 

— U.-L-L zc U.UL 






(-0.11 ±0.01) 


(-0.07 ±0.01) 


/dm - M dyn 




0.28 ± 0.02 


0.33 ±0.03 






(0.27 ±0.01) 


(0.26 ±0.03) 


M dyn /Mi, - 


Rett 


0.40 ±0.03 


0.16 ±0.04 






(0.38 ± 0.03) 


(0.12 ±0.05) 


M dyn /Mi, - 


n 


0.27 ±0.04 


-0.24 ±0.05 






(0.25 ± 0.05) 


(-0.13 ± 0.05) 


M dyn /Mi, - 


a 


0.78 ± 0.13 


0.76 ±0.14 






(0.69 ± 0.12) 


(0.61 ±0.13) 


M dyn /Mi, - 


Mi, 


0.19 ± 0.03 


0.10 ±0.02 






(0.13 ± 0.03) 


(0.04 ± 0.02) 


M dyn /Mi, - 


(P*) 


—0 17 + 01 


— 06 + 01 






(-0.16 ±0.01) 


(-0.05 ± 0.01) 


M dyn /Mi, - 


M dyn 


0.36 ± 0.03 


0.23 ±0.03 






(0.34 ± 0.02) 


(0.20 ±0.03) 


(PDM) — -Reff 


-1.91 ±0.06 


-2.16 ±0.09 






(-1.96 ±0.06) 


(-2.30 ±0.09) 


(pdm) — n 




-1.70 ±0.16 


-4.23 ±0.22 






(-1.77 ±0.14) 


(-3.25 ± 0.12) 


(pdm) - cr 




1.06 ±0.61 


1.88 ±0.74 






(0.77 ±0.43) 


(0.96 ±0.76) 


(pdm) - Mi, 




-0.91 ±0.11 


-0.78 ± 0.18 






(-1.08 ±0.08) 


(-1.28 ± 0.12) 


(Pdm) - (p*) 


0.70 ±0.02 


0.83 ±0.03 






(0.71 ± 0.02) 


(0.87 ±0.03) 


(pdm) — M dyn 


-1.05 ±0.08 


-0.68 ± 0.17 






(-1.12 ±0.08) 


(-0.93 ± 0.16) 



4 In each plot the sample is complete at log R c g/kpc ^ 0.5 kpc, 
log<r/(fcm/s) ^ 2.2, logM*/M > 10.5, log (p*)/(A/ Q /fcpc 3 ) sg 
8.5 and log M dyn / Mq ^ 10.4. Because of the large scatter in the 
diagram with Sersic n vs. SDSS Petrosian magnitude, we were 
not able to define a reliable completeness limit with respect to 
the n parameter. However, as shown in Fig. [3] (see also Fig. 0, 
the correlation of DM fractions (densities) with n arc approxi- 
mately linear over the entire range of this parameter, implying 
that the completeness limit is unimportant to characterize these 
correlations. 

5 We consider linear relations of the form Y = a + bX , where 
b is the slope. A bootstrap method is applied to estimate b and 
its uncertainty. In practice, we bin the data with respect to the 
variable X. For each bin, we randomly extract 50% of the points 
in that bin, computing the corresponding median value of Y. We 
perform 1000 iterations, each time computing the b from a linear 
fit of the median values of Y vs. the median values of X in all dif- 
ferent bins. The final value of b and its uncertainty are the median 
and <r of the distributions of slope values among all iterations. 



parameter space, studying /dm is not the same as studying 
M dyn /Mi,. In fact, considering the non-linear relation be- 
tween /dm and M dyn /Mt,, computing the average M dyn /Mi, 
and plugging it into the definition of /dm is not the same as 
computing directly the average /dm- 

In agreement with NRT10 and T+10, we find a tight 
and positive correlation between SIS /dm and Reg, which 
may be interpreted as a physical aperture effect, where a 
larger -R eff subtends a large r portion of a galaxy DM halo 
(see also lAuger et alj|2010bl ). A steep correlation also holds 
between DM fraction and Sersic n, namely, galaxies with 
steeper light profiles have higher central DM fractions. We 
also find that, independent of the adopted mass proxy (i.e. 
a, Mi,, Mdyn), more massive galaxies have the largest DM 
content (/dm ~ 0.6), while less massive ones have /dm~ 0.3. 

These findings are consistent with other results in 
the lite rature, based on different samples and method- 
ologies (|Padmanabhan et all |2004 ICappellari et al.l 120061 : 
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Figure 3. DM fraction within R c g as a function of A"-band i? e fj, uk, velocity dispersion cr, stellar mass M*, central average stellar 
density (/>*) and dynamical mass, M dyn . The symbols are as in Fig. [2] The vertical dotted lines set the completeness limit determined 
as discussed in the text. 



iHvde fc Bernardil l2009al : T+09; ICardone fc Tortoral l20ld : 
NRT10). A n exce p tion i s the conceptually related SDSS 
analysis by iGrilld l|201Ch , who found a fairly weak /dm- 
R E ft and /dm-o" correlations, and an anti-correlation between 
/dm and M*. These differences may be attributed to the 
use of a constant K, while our effective K values varied 
from galaxy to galaxy dependi ng on t heir in dividual Sersic 
parameters. We also note that iGrillol ([2010) discussed pro- 
jected rather than de-projected DM fractions. He converted 
the SDSS velocity dispersion to a R e g/8 circular aperture 
and adopted an expression for the dynamical mass similar 
to our Eq. [1] in his equation 4. The K = n factor in his 
equation 4 is equivalent to K = 2 if the dynamical mass 
is de-projected, and is slightly lower than the average value 
~ 2.15 we have reported above. There is also a difference in 
using 7f-band versus r-band effective radii (discussed be- 
low) and in other details related to the sample selection 
(such as the redshift range). Fig. [3] also shows a sharp anti- 
correlation between DM content and central average stellar 
density, which has not been reported in the literature so 
far. Galaxies with denser stellar cores have lower DM frac- 
tions (i.e. /dm ~ at p* ~ 3 x lO 9 M0kpc -3 ), while /dm 
values as high as ~ 0.8 are found at the lowest densities 
(p* ~ 10 6 M©kpc -3 ). This trend results from the fact that, 
on average, higher stellar densities correspond to smaller ef- 
fective radii, implying a lower /dm- For comparison with 



previous studies (see Sec. 12. 3p . Fig. also plots results ob- 
tained with the const-M/L model. The /dm values are still 
found to increase with R c s, but with shallower slope relative 
to the SIS, while the trend with the Sersic n is inverted. The 
trend with a is also similar to that obtained with the SIS, 
but at fixed a the /dm values are smaller by ~ 0.1-0.2. The 
const-M/L model implies no correlation of /dm with stellar 
mass, and a shallower trend (relative to the SIS) with both 
central stellar density and Mdyn- 

From the slope values reported in Table [JJ we notice 
that using either /dm or Md yn /M* as the dependent vari- 
able in the fits leads to different conclusions about the com- 
parison of SIS and const-M/L trends. For instance, the cor- 
relation of /dm with J? e ff is steeper for const-M/L than SIS, 
while the opposite holds for the Md yn /M*-R ff correlation. 
Whilst the use of /dm and Md yn /A/* is not truly equivalent, 
because of the scatter in DM content at each point of the 
parameter space (see above), the apparent discrepancy seen 
in Table [T] for, e.g., the /dm- and Mdy n /Af + -J? e ff relations, 
is indeed caused by the non-linear relation between /dm and 
Mdyn/M*: a difference in Mdyn/M* (such as that between 
SIS and const-M/L models) corresponds to a smaller dif- 
ference in /dm at higher, relative to lower, Md yn /A/* (i.e. 
going from low- to high-i? c g values). See also Fig. f3] where 
this effect is make clear by plotting /dm- and M dyn / M*-R e g 
relations. 
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0.8 1.0 1.2 
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Figure 4. M dyri / (top panel) and /dm (bottom panel) within 
R e ff in terms of i? e ff for a SIS and const-M/L. The plot is limited 
to the region where the sample is complete. The symbols are as 
in Fig. [2] and [3] 



The red curves in Fig. [3] also show how median /dm 
trends change when one excludes objects with Md yn < M*, 
i.e. /dm < 0. Except for the trend with n where a small 
rigid offset is seen (~ 0.05 for the SIS), all correlations are 
only affected at one end, i.e. at low mass, velocity dispersion, 
and R c ff, and at high (p*), where our sample turns out to 
be incomplete (see above) . Hence, the issue of negative /dm 
values only produce a slight effect on the slopes reported in 
Tab. [I 

Notice that negative /dm values can be caused by mea- 
surement errors on galaxy parameters (i.e. a, R B a, and M*), 
a failure of the mass model to estimate Md yn , and/or some 
systematics in the stellar mass estimates, such as the as- 
sumption of a given fMF (see Sec. I3.3fl . In order to assess if 
these effects can fully account for the negative /dm values in 
our sample, for each M*, we assign a mock Mdyn value, ac- 
cording to the Mdyn-Af* relation of Fig. [2] (for a SIS model). 
Shifting the M* and mock Md yn values according to the esti- 
mated statistical uncertainties on these quantities (~ 0.1 dex 
and 0.18 dex, respectively), we find ~ 2% of data-points to 
have Md yn <M*. This fraction rises to ~ 12% (i.e., what we 
actually measure for the SIS model), if stellar masses are 
overestimated by a factor of 1.6, which is within the factor 
of 2 systematic uncertainty on M* due to different effects 
(i.e. mainly the extinction law and IMF; see Sec. 12.21 and 
Paper V). 

An alternative explanation for the unphysical /dm val- 
ues is the choice of the galaxy mass profile, which could 
be inappropriate for (some) low mass galaxies. For instance, 
dynamical masses would become systematically larger (than 



p(r) oc r~ a 

with a < 2, is adopted (with M(r) oc r 3 ~ a , i.e. a mass 
profile shallower than the SIS). Different observations have 
argued that the (dynamical) mass profile of a galaxy might 
change as a function of its stellar mass. In particular, as 
already pointed out in Sec. 12.31 massive ETGs are quite 
well reproduced by nearly-isothermal profiles (e.g.lKochanekl 
Il99ll ; iKoopmans et al.ll2006l ; iGavazzi et ail 120071 ). while in 
some intermed iate-mass ETGs the mas s is found to fol- 
low the light ([Romanowskv et aL 20031; Napolitano et all 
120051 ; iDouglas et al.l 120071 ; iDeason et al l boi2[ ). Similarly, 
non-homology of the mass profiles has also been suggested 
by T+09. Qualitatively, this would make the trends of /dm 
in terms of Reff, a, M* and Mdyn steeper with respect to 
the case of a pure SIS model, and shallower (or inverted) for 
the cases of /dm versus stellar density and Sersic index. 

Whilst most of these results are qualitatively consistent 
with those found in the (recent) literature, one should no- 
tice that, in contrast to previous studies, our DM fractions 
are computed within the Tf-band (rather than optical) ef- 
fective radius. Hence, our /dm estimates are less affected by 
the existence of metallicity and age gradients in ETGs, dust 
extinction and (low fractions of) young stars. In order to 
permit a more direct comparison with results from our pre- 
vious work (T+09), we recomputed SIS /dm estimates using 
g-band (rather than K-h&nA) R e ff values. Fig. [5] compares 
the resulting gr-band median trends with the -B-band trends 
obtained from T+09. For the trend with stellar mass, we 
find good agreement between the present and our previous 
work. However, the plots of /dm vs. R c g and M dyn reveal 
some discrepancies. At low R c ff and A/dyn, the present sam- 
ple becomes incomplete and we miss galaxies with high /dm ■ 
This biases our trends towards lower /dm values with respect 
to T+09. At high J?cff, our i? e ff values are larger than those 
of T+09 (because of the different way structural parameters 
are estimated; see App.|X]), explaining the shift of the /dm 
vs. -Reff trend towards larger -R c ff with respect to T+09. On 
the other hand, the present /dm vs. Mdyn trend is shifted 
towards lower Md yn values relative to T+09, implying that 
some further systematic differences exist between the Md yn 
estimates of the two datasets. Understanding the origin of 
such differences is far from trivial as we do not have galax- 
ies in common between the two samples. Differences might 
exist because of different sample selection criteria, and/or 
different methods used to estimate velocity dispersions. 



3.3 Constraints on IMF 

The present analysis is based on the assumption that the 
stellar IMF is universal, at least within the parameter range 
covered by our sample. In principle, for each given galaxy, 
one could adjust the IMF so that M*=Md yn (or, equiva- 
lently, /dm= 0). However, since /dm correlates with various 
galaxy parameters, any adjustment of the IMF would imply 
that the IMF itself changes syst ematically with these p aram- 
eters, in the same way as /dm (|Cappellari et al.ll2012l ). Mo- 
tivated by the recen t claim that the IMF depends on galaxy 
velocity dispersion (jvan Dokkum fc Conroy|[20Tll . hereafter 
vDCll), with high- (relative to low-) a ellipticals having 
a bottom-heavier IMF, we explore to what extent one can 
change the IMF in order to fulfill the condition /dm= 0. 
We estimate the expected variation of stellar mass, 
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log R ef f [kpc] log M 4 [M Q ] log M dyn (R eff ) [M ] 



Figure 5. /dm within g-band R c ff in terms of R e ff, M* and Md yn for a SIS. Solid curves and shaded regions are for our sample, as in 
Fig. [2] while points with error bars are from T+09 and NRT10. The vertical dotted lines mark the completeness limits as in Fig. [3] 
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Figure 6. Trends of M^ yn /M ir as a function of X-band R c g, stellar mass M* (assuming a Chabrier IMF), velocity dispersion, and 
dynamical mass, M,j yn . Black curves and grey regions are the same as in Figs. [2] and [3] Horizontal lines correspond to the relative 
variation of stellar mass, M* imf/-M* - with respect to a Chabrier IMF - when adopting different IMFs, with slopes 1.35 (i.e. a Salpeter 
IMF; blue), 1.85 (green), and 2.05 (cyan). The red line corresponds to the case of a Chabrier IMF (M^imf = M+). The vertical dotted 
lines set the completeness limit as in Figs. [3] and [5] For /dm= mp/M t , one can find the IMF slope that would make M* equal to 
Md yn , resulting in a galaxy with no DM. 



M+jmf /Af,, relative to a Chabrier IMF, when different 
IMFs are adopted. Here, M* is the stellar mass estimated 
with a Chabrier IMF, while M*,imf is the stellar mass we 
would estimate with a different IMF. In practice, we con- 
sider three power-law IMFs, with slopes 1.35 (i.e. Salpeter), 
1.85, and 2.05 (i.e. a very bottom-heavy IMF, as suggested 
by vDCll for high-cr ellipticals). The M*,imf/M* is esti- 
mated as the ratio of the stellar M/L between two SSPs 
having a power-law and a Chabrier IMFs, respectively. To 
compute the stellar M/Z0, we adopt the BC03 synthesis 
code, for old (10 Gyr) SSPs, with solar metallicity. 

Fig. [6] plots the Md yn /A/* trends as a function of R c ft, 
M+, Mdyn, and a (black curves), with horizontal lines mark- 
ing M^imf/A/* for the different IMFs. The intersections of 
the horizontal lines with the black curves define the values 
of R c a, M+, Mdyn, and a, for which a given IMF slope would 
imply /dm= 0. The Figure shows that in order to account 
for the apparent trend of /dm with _R e ff (and Md yn ), galax- 
ies with the largest radii (and dynamical masses) should 
have an IMF slope as steep as 2.05, while at the lowest 



6 We compute the M/L values in the K"-band, as this is less 
sensitive to the presence of young stars and metallicity of a stellar 
population. 



R e fi (and Md yn ) a Chabrier (or even a bottom-lighter) IMF 
would be required. Interestingly, at high velocity disper- 
sion (loga ~ 2.45), we see that the IMF slope cannot be 
larger than ~ 1.8, which seems to contrast with the finding 
of vD Cll, but is in good agreement with ICappellari et al.l 
|2012l ) (see their no-DM case). 

The exercise described so far shows that the estimated 
trends of /dm might be significantly biased if the IMF 
varies substantially as a function of different galaxy pa- 
rameters. However, it remains unclear if and (more impor- 
tantly) why such correlated variation of IMF slope with dif- 
ferent parameters should be present. A nother key to the /dm 



trends might be a l ternative gravity JCardone et al. 2011 
Lubini et all l201ll ; iRichtler et~aH 1201 ll ; iNapolitano et al 



al] 

20121 ). which, anyway, is beyond the scope of the present 
work. 



3.4 Dark matter densities 

Following T+09 and NRT10, we also analyze the average 
central DM density of ETGs, (pdm)- The (pom) values 
may be considered the fossil record of the ambient density 
at the time of initial halo collapse, modulo possible mass- 
redistributing interactions between DM and baryons. While 
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Table 2. Comparison of the /dm — X and /dm — (Y) relations 
(black and red curves in Fig. [8] see the text). We flag with 'Yes' 
('No') those cases where /dm — (V) overlaps (does not overlap) 
with /dm - X. 







RaS 


Y> 


n 


a 


<P*> 


X = 


log R eS 




No 


No 


No 


Yes 


X = 


logM* 


No 




No 


Yes 


Yes 


X = 


n 


Yes 


No 




No 


Yes 


X = 


logo- 


No 


No 


No 




No 


X = 


log (fii,) 


Yes 


No 


No 


No 





/dm is defined relative to the baryon mass, (pdm) quanti- 
fies the properties of the DM component alone, allowing for 
more direct comparisons to cosmological models (see Sec.0. 

Fig. [7]is the same as Fig. [3] but plotting (pdm), rather 
than /dm, as a function of structural parameters and galaxy 
mass proxies. The correlations are again fitted with lin- 
ear relations, taking into account the completeness limit 
of each diagram, with the slopes reported in Table [T] We 
find that (pdm) decreases with Reg, n, M+, and Md yn - The 
best-fittin g power laws are con sistent with res ults from the 
literat ure ((Thomas et afll2009t T+09; NRT10; [Auger et all 
l2010bl ; T+10). On the other hand, (pdm) 

increases with a 

and (p*). All trends are qualitatively independent of the 
adopted mass model, although in general the const-M/L pro- 
file produces steeper correlations than the SIS. As shown in 
Fig. at low R cff , n, a, M dyn , and high (p*), (pdm) is 
model independent, while the discrepancy between SIS and 
const-M/L models is larger at high R c fi, n, a, Md yn , and 
low (p*), making the const-M/L correlations steeper. This 
is due to the increasingly important role of mass profile ex- 
trapolation in these regimes. 

Under the assumption that DM density profiles are uni- 
versal, at least for all galaxies in our samples, we confirm the 
result of NRT10, that DM profiles in ETGs are very cuspy, 
with pdm{t) ~ r~ 2 . This conclusion is reached as follows. 
For a power-law DM density distribution, pDAi(r) oc r~ a , 
we have Mom = r' i ~ a (for a < 3 ). It follows that 
(PdmH-Rch) oc R cH ~ 3 M D M(R e ff) oc R ca ~ a . Since the slope 
of the (pdm) - Rcff correlation is ~ 2 for both SIS and 
const-M/L, we conclude that ot ~ 2, independent of the 
adopted mass model. This conclusion relies on the assump- 
tion that trends in the average values of (pdm) reflect trends 
in th e mass profiles of individual galaxies (cf. IWalker et al.l 
120091 ; NRT10). We will return to further implications of this 
assumption in Section [7] 



4 "DARK MATTER PLANE" OF ETGS 

Some of the correlations discussed so far, between DM frac- 
tion (density) and different galaxy parameters, could in prin- 
ciple be secondary in nature, arising simply from mutual 
correlations among the galaxy parameters themselves. For 
instance, given the fuM-R e ff relation, a correlation between 
/dm and stellar mass is also expect ed to exist, just b ecause of 
the size-stellar mass relation (e.g. lShen et al.ll2003l ). There- 
fore we seek to determine the fundamental parameters that 
relate to the DM matter content of ETGs. We follow the 



same procedure adopted in lLa Barbera et all (|2010dl ) (Pa- 
per IV) to establish the main drivers of correlations among 
internal colour gradients and other observed parameters of 
ETGs. As illustrated in Fig.[8l we consider a given galaxy pa- 
rameter, X, out of the quantities {R c ff, n, a, M*, and (p*)} 
(from top to bottom in the Figure), and the corresponding 
correlation with /dm, /dm(^) (black curves in each panel). 
For a given bin of X, we compute the median value ((Y)) of 
another galaxy parameter, Y , with F/I, and then, given 
the /dm (Y^j correlation, we compute /dm((^}) (red curves 
in the Figure). The /dm({^ / )) is the correlation between /dm 
and X expected from the correlation between /dm and Y, 
and the correlation between X and Y themselves. Each row 
in Fig. [8] from top to bottom, corresponds to a given X, 
while the panels from left to right show the /dm((V)) (see 
red curves) by varying Y . If a given /dm((V)) differs from 
/dm(X) (i.e. the black and red curves in the corresponding 
panel do not overlap), we can conclude that the correlation 
between /dm and X is more than just a reflection of one be- 
tween /dm and Y. We evaluate by visual inspection whether 
or not two curves differ, as summarized in Table O for each 
pair of X and Y. We find that the correlation between /dm 
and X — log Ah (panels 2a-2d) is equivalent to that ex- 
pected from either /dm and a (panel 2c), or /dm and (p*) 
(panel 2d), i.e. stellar mass is likely not a genuine driver of 
the DM fraction in ETGs. In the same way, the correlation 
with Sersic n is equivalent to that expected between /dm and 
Retf (or log(p*)) (panels 3a and 3d). On the contrary, the 
correlation between /dm and a is a genuine one, as all the red 
curves in panels 4a-d of Fig. [S] differ, being flatter than the 
observed /DM(logcr). Finally, we see that the /DM(log -R c fr) 
and /dm (log (p*)) correlations are completely equivalent to 
each other [j, and are not caused by any correlation of /dm 
with some other parameter. Since p* is computed from Rett, 
we conclude that the main parameters driving DM fractions 
are galaxy size and central velocity dispersion. The same 
result holds when repeating the above analysis using DM 
density, rather than /dm- 

Since there are two observed parameters driving the 
DM content of ETGs (R B g and a), we analyze here the cor- 
relation of DM with both quantities, using SIS-based DM 
estimates. Here, we adopt the logarithmic ratio of dynami- 
cal to stellar mass, log (Md yn /A/*)- This avoids the issue of 
negative /dm values and can provide a more direct connec- 
tion to the FP relation (see § [5}. We consider a DM plane 
relation: 

log A ^ yn = cidm log Rcff + boM log <r e + com, (2) 

where cidm and 6dm are the slopes, and cdm is the offset. 
We perform a linear best-fit of log (A/d yn /A/*) vs. log_R c ff 
and a(R e ff), minimizing the sum of absolute residuals along 
the M dyn /M+ axis. This gives a D M = 0.380 ± 0.008, 6dm = 
0.647 ± 0.03, cdm = —1.38 ± 0.06 (where errors are la sta- 
tistical uncertainties). Both «dm and 6dm are significantly 
different from zero, implying that DM content is indeed a 
function of both R e ff and a, consistent with what we found 



7 In fact, the observed /DM(logfl c ff) relation is fully consistent 
with that expected from /DM(l°g (p*)) (panel Id), and vice- versa 
(panel 5a). 
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Figure 7. Average central DM density as a function of _R"-band il e g, ra, velocity dispersion cr, stellar mass M*, central average stellar 
density (p*) and dynamical mass, M dyn . The symbols are as in Figs. [2] and [3] except as for the red curves in that plot, here only the 
galaxies with M dyn > M* are used. The vertical dotted lines set the completeness limit as in Figs. [3] and [5] 



above (i.e. that the correlation of /dm and 7? e ff is not equiva- 
lent to that with a, and vice- versa). Fig.[9]compares the (log- 
arithmic) correlations of Md yn / vs. R B a (top), Md yn /M+ 
vs. a(R c ff) (middle), and those between Md yn /M+ and both 
7? e ff and a(R e ff) (i.e. the DM plane; see bottom panel). See 
the slopes reported in Table [1] For all panels, the red lines 
are the best-fitting relations obtained by the same kind of 
fitting procedure, i.e. minimizing the sum of absolute resid- 
uals along the Md yn /M* axis. A robust estimate of the rms 
along the Md yn / 'M* axis is also reported in the Figure for 
each correlation. Notice that the scatter of the DM plane 
(rms = 0.149 ± 0.002) is smaller (albeit by a few percent) 
than that of the M dyn /AL - a(R cB ) (rms = 0.200 ± 0.002) 
and Mdyn/M* - R c a (rms = 0.165 ± 0.002) correlations. 

Projecting the Md yn /M+ rms along the log J? e fF axis, we 
find a scatter of ~ 0.4 dex, i.e. much larger than that of the 
FP relation (~ 0.1 dex; see Paper II). This large dispersion 
implies that the coefficients of the DM plane depend signif- 
icantly on the fitting procedure one adopts to derive them. 
For instance, minizing the rms of absolute residuals along 
logi? e ff, we find «dm ~ 0.84 and 6dm ~ 0.78, which are sig- 
nificantly different from those obtained from the Md yn / 'M+ 
best-fit (see above). The scatter of the DM plane decreases 
to ~ 0.25 dex along logi? c ff, still significantly larger with 
respect to the FP. 



5 DARK MATTER FRACTIONS AND 
FUNDAMENTAL PLANE 

The correlations of DM content with galaxy parameters 
might help to shed light on the origin of scaling relations 
of ETGs, i.e. the tilt of the FP relation. The FP can be 
written as 

rogi? c ff = a_Fplogcr + &Fp(A t *)e + constant, (3) 

where cifp and 6fp are the slopes, and (^t*) e is the mean sur- 
face brightness within R c s . Under the assumption of homol- 
ogy, the tilt can be parameterized as a variat ion of Md vn /L 
with luminosity or mass, i.e. Ma yn /L oc L a (|Dressler et al.l 
Il987t ). or M dyn /L oc M dyn 7 (with 7 = a/(l + a)). The M/L 
can be rewritten as 

Afdyn Mdyn 

L K L ' { > 

where M+/L is the stellar mass-to-light ratio. In the K- 
band, stellar p opulation effects make a negligible contribu- 
tion to the tilt (|La Barbera et al"1l2010bh , implying that non- 
homology and/or a variation of DM content (i.e. Mdyn/M* 
changing with luminosity /mass) should explain the tilt. 
Fig. [3] shows that the central DM content of ETGs does 
actually change significantly with mass and luminosity. In 
particular, we find M dyn /M* oc M dyn °' 36±0 - 03 , or 7 = 
0.36 ±0.03 (see Table[l|. This value is somewhat larger than 



12 Tortora et al. 



0.8 
0.6 



tr 
2 0.4 



0.2 
0.0 



0.8 
0.6 



g 0.4 



• Qa) ^^^=^^ ■ 


■ ^^^^^^^^ 


. (1c^ ^^j^ . 









0.0 0.5 1.0 1.5 
X = log R elf [kpc] 


0.0 0.5 1.0 1.5 

X = log R eff [kpc] 


0.0 0.5 1.0 1.5 

X = log R ef f [kpc] 


0.0 0.5 1.0 1.5 

X = log R ef f [kpc] 


■ (2a) Y = log R efl ■ 


. (2b) Y = n . 


■ (2c) Y = tr . 


■ (2d) Y - log <p*>- 


10.0 10.5 11.0 11.5 

X = log M, [M Q ] 


10.0 10.5 11.0 11.5 

X = log M, [M Q ] 


10.0 10.5 11.0 11.5 

X = log M, [M ] 


10.0 10.5 11.0 11.5 

X =log M, [M e ] 


■ (3a) Y = log R eff ■ 


■ (3b) Y = log M, . 


■ (3c) Y = tr . 


■ (3d) Y = log <p,>- 


















4 6 8 10 

X = n 


4 6 8 10 

X = n 


4 6 8 10 

X = n 


4 6 8 10 

X = n 



0.2 
0.0 



0.8 
0.6 



DC 

= 0.4 



0.2 
0.0 



0.8 
0.6 



1 0.4 



. (4a) Y = log R e „ ■ 


• (4b) Y = log M, . 


■ (4c) Y = n 


■ (4d) Y = log </>,>■ 






2.0 2.2 2.4 

X = log tr 


2.0 2.2 2.4 

X = log tr 


2.0 2.2 2.4 

X = log tr 


2.0 2.2 2.4 

X = log tr 



0.2 
0.0 



0.8 



■ (5ak 


Y = log R e „ • 


■ (5b) 


Y = k 


>g M, ■ 


• (5c) 


Y = n 


■ (5d) 


Y = logtr • 





















0.6 












0.4 






0.2 



0.0 



4 6 8 10 4 6 8 10 4 6 8 10 4 6 8 10 

X = log <p,(R efl )> [M Q /kpc 3 ] X = log <p,(R eff )> [M Q /kpc 3 ] X = log <p»(R efr )> |M Q /kpc 3 ] X = log <p,(R eff )> [M Q /kpc 3 ] 



Figure 8. Tests performed to establish the observed ETGs' parameters driving the correlations with central DM fraction. Each panel 
shows the observed correlation between /dm and a given parameter X (black solid curve), as well as the correlation expected from that 
of /dm and another galaxy parameter Y (i.e. Y ^ X), and that between Y and X (red solid curve). Each row in the plot corresponds to 
a given X, out of Reff, M+, n, tr, and (p*) (from top to bottom). For a given row, the panels show the same fnwi-X correlation (black) 
and that expected by varying Y (from left to right). For each panel, the corresponding Y is reported in the upper-right corner. 



that found by some previous studies, i . e. 7 ~ 0.15-0.20, 
(|Hvde fc Bernardil 120093 . iGallazzi et "aD 120061 and T+09), 
possibly due to differences in the way structural parameters 
are estimated from different sources (see App.fX]). The value 
of 7, as obtained by fitting the Md yn /M + vs. Md yn relation, 



is close to that inferred from FP coefficients, using the virial 
theorem un der the assumption of h omology (7 ~ 0.2; see, 
e.g., T+09; lLa Barbera et alll2010bh . This means that DM 
content variation might be able by itself to explain the en- 
tire fraction of the tilt which is not due to stellar population 
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Figure 9. DM fundamental plane. Top Panel, log M^ yn / M* as 
a function of R e g. Middle Panel, log M^yn/M* as a function of 
cr e . Bottom Panel, log M^ yn /M i , as a function of the combination 
iDM log ii c ff + 6dm log cr e + cdm- The cyan lines are the medians, 
while the red ones arc the best fits. 



effects (i.e. a variation of age and/or metallicity with mass). 
We also notice that, in the context of the const-M/L model, 
the variation of DM content with mass is flatter than with 
the SIS model, which would lead one to conclude that both 
non- homology and a change of DM content contribute to the 
tilt (see T+09 for details). Last, but not least, one should 
notice that connecting the slope of the Mdyn/Af* vs. Md yn 
relation to the tilt of the FP is not trivial. The dynamical 
mass entering the virial theorem is actually the total mass of 
a galaxy, i.e. all the gravitationally bound matter of the sys- 
tem. Hence, the difference of FP coefficients from the virial 
theorem expectation might be more related to the total, 
rather than central (< 1 R c a), DM content of ETGs. Unfor- 
tunately, estimating total, rather than central, DM fractions 
requires a large extrapolation of the central velocity disper- 
sion, much more dependent on the mass model one assumes 
to infer our default Md yn values. With this caveat in mind, 
we try here to establish a connection between the DM plane 
and the FP, combining Eq. ([2]) with the following expression 
of the virial theorem, obtained under the assumption of ho- 
mology (whereby kinetic energy is oc Ma yn a 2 and potential 
energy is oc Mj yn /R cf i): 



where G is the gravitational constant. We obtain an equation 
formally identical to the FP: 



logi? c ff = a log <7 + b(E*) e + constant, 



(6) 



with a = (2 — 6da/)/(1 + om/), b = 0.4/(1 + <2dm), and 
(£*} e is the mean stellar mass density within R c g. 

Using the values of aoM and bom from Sec. [4] we ob- 
tain a = 1.00±0.02, b = 0.290=10. 002. These coeffici e nts can 
be compared to those derived by iHvde fc Bernardil l|2009bh 
for the stellar mass FP, i.e. a* ~ 1.4 and 6* ~ 0.32 (when 
using a direct fitting method, see their tab. 3), and those 
we derived for the K-b&nd FP in Paper II, i.e. ax ~ 1.55 
and bK ~ 0.32. While 6 is very consistent with 6* and bx , 
a is smaller than both a* and clk- This discrepancy might 
be explained by several effects, such as selection effects (see 
e.g. Paper II), fitting procedure (see comments about DM 
plane fitting procedure in Sec. |4j, and the method (i.e. the 
model) used to estimate DM fractions. In general, the fact 
that a and b, as derived above, are significantly smaller than 
the homologous expectations of the virial theorem (i.e. a = 2 
and b = 0.4), and similar (or even smaller) with respect to 
the observed coefficients of the stellar mass (or A'-band) 
FP, may lead one to conclude that most of the tilt is due, 
indeed, to a variation of DM content with galaxy mass (or 
cr and Red, as shown here). However, we still remark, as 
noticed above, that the DM plane is a correlation among 
central quantities of ETGs (with dynamical mass being es- 
timated within 1 Rch), while the virial theorem is a global 
relation (between kinetic and potential energy within an in- 
finite aperture). Thus, the above conclusion about the origin 
of the tilt remains uncertain. On the other hand, the exis- 
tence of a central-DM plane is a robust result (see previous 
section), and may be explained by the DM halos of ETGs 
being only approximately universal. A larger R e g encloses 
a larger portion of the halo, implying /dm to increase with 
R e g (as discussed by, e.g., NRT10). At given R c g, a larger 
cr means a deeper central potential well, i.e. more DM in 
the centre, implying that Mdyn/ 'M* increases with both R c s 
and a. 



6 CENTRAL DARK MATTER AND GALAXY 
ENVIRONMENT 

In a hierarchical paradigm of galaxy formation, the (dark) 
matter content of the outer, less bound regions of a galaxy 
can be stripped off, as the galaxy is accreted into a bigger 
halo, becoming a satellite. If stripping is effective enough, 
one might see some variation also in the central DM mat- 
ter fraction of galaxies, as a function of the environment 
wherein they reside. Environment-driven interactions, such 
as (both major and minor) merging with other group mem- 
bers, may further change the DM fractions of galaxies re- 
siding in groups. In this section we investigate the effect of 
galaxy environment on the central DM fractions of ETGs. 
Fig. [10] plots /dm as a function of both structural parame- 
ters and mass proxies (as in Fig. [3J), for ETGs classified as 
field and group galaxies, the latter being splitted between 
satellites and centrals (see Sec. [2|. We show here only SIS 
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Figure 10. Environmental dependence of the correlations between DM fraction and structural parameters/mass proxies. Only the SIS 
model is shown. Solid black lines and shaded regions are the same as in Fig. [2] and refer to the entire sample. Blue, red, and green lines 
show the median trends for field, satellite, and central galaxies, respectively. 



results, as const-M/L trends would not add any relevant in- 
formation to the analysis. 

We find that field and satellite galaxies exhibit fully 
consistent /dm trends while somewhat different trends are 
detected for centrals. In fact, when compared to the average 
trends of field and satellite ETGs, centrals have, at fixed 
n (Mdyn), larger (lower) SIS /dm values; while all remain- 
ing correlations do not show any significant environmental 
dependence (although at low a centrals also tend to have 
lower /dm). These differences are just due to the fact that 
for centrals the M* (f? c ff) is typically larger relative to field 
and satellite galaxies, as seen by the M* and i? e fi ranges for 
central galaxies (green curves) in the /dm~M* and fuM-R e B 
diagrams of Fig. [TO] Since /dm increases with both M* and 
R e H, centrals tend to have, on average, larger /dm, while at 
fixed Mdyn, a larger M* implies a lower /dm- In order to ac- 
count for a possibly different shape of the DM halo of satel- 
lite (relative to field) ETGs, we have also computed their 
/dm values by using a truncated SIS profile (see Sec. 12.31) . 
We find no significant difference with respect to the results 
of a pure SIS model. We also test for environmental depen- 
dence in the DM plane of ETGs. Minimizing the rms of abso- 
lute residuals along Mdyn/M*, as done in Sec.[4]for the entire 
sample, we obtain odm = 0.39±0.02 (&dm = 0.58±0.05) and 
a D M = 0.35 ± 0.01 (&dm = 0.60 ± 0.04) for field and group 
ETGs, respectively, i.e. we still do not detect any significant 
environmental dependence. In summary, we do not find any 
significant environmental variation of the correlations be- 
tween central DM content and structural parameters and 
mass of ETGs. This result can be discussed in light of what 
we found in Paper III, i.e. that the tilt of the FP is larger 
for group, relative to field, ETGs. Since we found this effect 



to be independent of the waveband where galaxy structural 
parameters are measured (from g through K), we concluded 
that it is explained by some "wavelength-independent" ef- 
fect, like a different variation of DM content with mass for 
group (relative to field) ETGs. This latter conclusion seems 
to contrast with what we have found here (Fig. llO|) . However, 
as pointed out in Sec. 21 in the present study we analyze cen- 
tral DM fractions, while the tilt of the FP may be more re- 
lated to the behavior of global quantities. Alternatively, the 
environmental variation of the tilt should be explained by 
a change of dynamical structure (i.e. non-homology) rather 
than DM content with environment. 



7 COSMOLOGICAL MODELS AND 
SIMULATIONS 

In order to interpret the observational results, we follow 
the same approach as in T+09, NRT10 and T+10, con- 
structing a set of toy-galaxy mass models, whose DM den- 
sity profiles are based on ACDM cos mological simulations, 
i.e. a iNavarro. Frenk fc White! |l99l hereafter NFW) pro- 
file. A different set of models is also constructed by apply- 
ing a suitable recipe for adiabatic contraction (AC) f rom 
baryon settling to the NFW profiles (|Gnedin et al. 1120041 ). In 
both cases (i.e. NFW and AC-NFW models), a de-projected 
Sersic law is used to describe the density profile of the stellar 
component. 

The NFW halos are ass umed to follow an a verag e mass- 
concentration relation, as in lNapolitano et all (|2005h . while 
the Sersic profiles are assumed to follow the _R e ft-M* and 
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Figure 11. Comparison of the DM content of ETGs to the predictions of galaxy-toy models, /dm an d (pdm) are plotted as a function 
of R c ff , Mi, and M<j yn in the upper and lower panels, respectively. Black and red curves, as well as grey shaded regions are the same as 
in previous plots. The solid and dot-dashed green lines are toy- model predictions for a AC-NFW and NFW DM halo (see the text). 



n — M* relations from our datc[f|. The virial DM mass 
is parameterized in terms of a star format ion efficiency 
esF = M*/(fi ba fM t ot), where f2 bar = 0.17 (|Spergel et all 
12007ft is the baryon density parameter and Mtot is the to- 
tal mass. We assume that esF v aries with M*, following the 
recipes of lConrov fc Wechslerl (|2009D . where esp decreases 
from a maximum of ~ 0.19 at log M*~ 10.4 down to ~ 0.004 
at logM*~ 11.8. These assumptions allow us to parametrize 
all toy-model properties (e.g. the central DM fraction) as a 
function of one single quantity (e.g. i? e ff, M+, or Md V n)- In 
Fig. 111! we compare the observed correlations of DM frac- 
tion (upper panels) and density (lower panels) with 7? e ff, 
Mi,, Mdyn (from left to right in the Figure), to the toy- 
models' predictions, where solid and dashed green curves 
correspond to AC-NFW and NFW models, respectively. 
Despite some differences in the trends, we find qualitative 
agreement between the SIS-based observational results and 
the toy-models, as the latter occupy a region similar to that 
of the data in each diagram. On the other hand, adopting 
an extreme const-M/L model would produce disagreement 
between data and theory, which predicts more DM on scales 
of ~ 10 kpc and beyond. 

8 The R a ff — Mi, relation is well reproduced by the best fit- 
ted relation logi? c ff = -8.7 + 0.86 log M,/M Q , while is logn = 
-5.3 + 0.58 log M*/M and logn = 0.27 + 0.05 logA/*/M for 
galaxies with log Mi, / Mq < 10.5 and J>10.5. 
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Figure 12. Comparison of the /dm~M 4 relation to the predic- 
tions of cosmological simulations. The black and red curves, as 
well as the grey shaded region, are the same as in previous Fig- 
ures. Red and blue dots correspond to simul ated brightest cluster 
galaxies from iRuszkowski fc Springell |200Sh . for contracted and 
uncontracted models, respectively. Black dots and open squares 
are th e results of hydrodynamical simulations of lOhorbe et al.l 
||2007|) , adopting high and low gas-to-stars conversion efficiency, 
respectively. 
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In more detail, uncontracted (i.e. NFW) models predict 
on average smaller /dm than contracted models (NRT10; 
T+10). Both sets of toy-models give similar /dm and (pdm) 
predictions at high mass (R c s), while they deviate more 
from each other at lower mass (R e a), with differences up 
to A/dm ~ 0.2 and A log (pdm) ~ 0.5 Mq kpc -3 . 

The AC model matches the observations well over a 
broad range of galaxy sizes and masses, which is remark- 
able since there was no fitting or fine-tuning involved. The 
NFW model fares less well, particularly for galaxies with 
small masses and sizes. As discussed in Section \3. 41 the ap- 
parent implication is that the central DM profiles have den- 
sity cusps that are steeper than in an NFW profile, and are 
broadly consistent with adiabatically-contracted NFW pro- 
files (see also NRTIO; T+10). Such conclusions will depend 
on the modelling ingredients, such as the IMF adopted, the 
esF vs Mi, relation etc., which it is outside the scope of this 
paper to examine in detail. Here we can only remark that 
one way to fix these discrep ancies is to allow the 6 sf to grow 
with mass in contrast to the lConrov fc Wechsler] (|2009l ) pre- 
dictions. 

In Fig. [TJ] we compare our /dm-M 4 relation with the 
results of two sets of cosmological simulations: (I) simu- 
lated brightest cluster galaxies from lRuszkowski fc Springell 
(|2009l ). using uncontracted and contracted DM models (blue 
and red dot s, respectively ) ; and (2) hydrodynamical simu- 
lations from lOnorbe et al l (|2007l . hereafter ODS07), for the 
cases where the authors adopted either a high (black dots) or 
low (black empty squares) gas-to-stars conversion efficiency. 

We find that contracted models are in excellent agree- 
ment with the observed SIS-based /dm values, while uncon- 
tracted models produce results closer to const-M/L /dm esti- 
mates, which is qualitatively consistent with our toy-model 
results. The SIS /dm~M* relation is also remarkably con- 
sistent with the simulations of ODS07, in the case of high 
conversion efficiency. On the other hand, const-M/L /dm's 
are more consistent to the predictions of models with low 
conversion efficiency. This is due to the fact that the low 
efficiency simulations of ODS07 are those producing lower 
/dm's. In fact, as noticed by the authors (see their Sec. 6.2), 
adopting a lower conversion efficiency yields lower _R e ff val- 
ues, which enclose a smaller portion of a galaxy DM halo, 
implying lower /dm values. Finally, one should notice that 
whilst our average /dm~M* relation is fully consistent with 
that of ODS07, the measured scatter around it is signifi- 
cantly larger than that of the simulations. This might be due 
to measurement errors, as well as the effect of some other pa- 
rameter which is not taken into account by simulations (e.g. 
galaxy stellar populations). Moreover, consistently with the 
results in Fig. III! it is encouraging that our 7? e ff-M* relation 
is broadly consistent with the relation predicted by ODS07 
for high efficiency models. 



8 CONCLUSIONS 

We have discussed the central DM content of massive, 
nearby (z < 0.1) ETGs, using high quality data from the 
SPIDER survey, including stellar masses estimated from 
optical- to NIR-photometry. Using the Jeans equations, we 
model each galaxy to estimate its dynamical mass, Md yn , 
within I i? fi in the K-ha,nA. We have adopted a SIS and, 



as a comparison, a non-homologous const-M/L model for 
the total mass profile, in order to investigate the system- 
atics induced by different assumptions for the mass profile. 
The recovered DM fraction, /dm, or equivalently the ratio 
Mdyn/A/* and the average densities, (pdm}, have been an- 
alyzed as a function of structural parameters and different 
mass proxies, comparing results with predictions of both toy- 
models and cosmological simulations. The main results are 
the following: 

• For a Chabrier IMF we find that only small fractions 
(~ 12 and 24% for the SIS and const-M/L, respectively) of 
galaxies have negative DM fractions, while the fractions in- 
crease to ~ 55 and 78% when a Salpeter IMF is adopted. 
Thus, under the assumptions that the total mass profiles 
are well described by a SIS or a const-M/L and the IMF is 
the same for all galaxies, our data favour a Chabrier (or a 
Kroupal200"ll ) IMF, in agreement with, e.g.. lCappellari et al.l 



20061 . On the other hand, relaxing the assumption of a uni- 



versal IMF, for the most massive galaxies in our sample 
(with mass comparable to the most massive early-type grav- 
itational lenses in the SLAC S survey) a Salpeter IMF is 
consis tent with observations (|Treu et alj|2010t lAuger et al.l 
l2010al ). We remark that, at low mass, a Chabrier IMF 
is al s o preferred for lat e-type galaxies (|Sonnenfeld et all 
l201ll ; iBrewer et al]|2012l ). and a Salpeter IMF is ruled out 
for ETGs, compared to a Ch abrier IMF, based on lens- 
ing and stellar popu lation data (|Ferreras. Saha. Burlesll2008l : 
iFerreras et~aHl2010h . 

• DM fractions and densities are, on average, higher when 
an SIS is adopted, when compared with the results using 
const-M/L profile. This result is also confirmed by the de- 
terminations of the virial coefficient in Eq. (JTJ, K, where, 
on average, K S is > K M/h . 

• The DM fraction is a steeply increasing function of 
R c g , Sersic i ndex, a, M+ and Md vn , c onfirming previ- 



ous fi ndings ([padmanabhan et"al] 12004 ; Cappellari et al.l 
20061; iHvde fc Bernardil l2009al; T+09; lAuger et all 1201 Obi; 



Cardone fc Tort ora 2010). Galaxies with denser stellar cores 



have also lower /dm values. The assumption of a non- 
homologous const-M/L profile yields weaker trends. The 
mass dependencies may provide important clues to the FP 
tilt. Assuming that central DM trends can be combined with 
the virial theorem equation, the SIS results point to a dom- 
inant role of DM. If a non homologous const-M/L profile is 
assumed, then almost all the tilt would be due to a non- 
homology. We have also verified that these results are quali- 
tatively consistent with our previous determinations (T+09 
and NRT10). 

• Similarly, the central average DM densities are decreas- 
ing functions of -R e ff, n and masses. Steeper trends are found 
when a const-M/L profile is adopted. In particular, (pdm) 
R c g ~ 1,8 and R c g ~ 2,2 for the SIS and const-M/L respectively, 
due to the steeper DM density profile for const-M/L at R c f(. 

• We have compared our results with predictions from 
ACDM toy-models, finding good qualitative agreement. As 
in NRT10, we have found that when a Chabrier (or Kroupa) 
IMF is adopted, a contracted NFW is preferred, while a 
bottom-heavy IMF like Salpeter would be more consistent 
with an unco ntracted NFW ( see, e.g. ICardone fc Tortoral 
|2010| ; NRT10; iTreu et alj|20ld ). 

• The role of the environment has been investigated after 
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classifying the galaxies in the field and in groups, and in 
the latter, centrals have also been examined explicitly. We 
find no difference in the estimated /dm between field and 
group galaxies, and only slight differences between centrals 
and non-centrals. This result is possibly due to an increased 
role of mergers for the central galaxies. 

We have finally shown that a very tight relationship 
similar to the FP exists between DM fraction, velocity dis- 
persion and size. We have determined the best fitted pa- 
rameters in Eq. @ and have inferred that DM might be the 
dominant driver of FP tilt, at least when a SIS is adopted, 
and under the assumption that central DM trends can be 
combined with the virial theorem expectation. In a forth- 
coming paper we will discuss the results introduced here 
in terms of the stellar population parameters derived from 
the fitting of spectra, by means of different synthetic pre- 
scriptions. In particular, within the general framework intro- 
duced in NRT10 and T+10, we will investigate correlations 
with stellar populations parameters such as age, metallic- 
ity, and alpha-enhancement. With further comparisons to 
simulations, these correlations may be used to constrain the 
evolution of stars and DM in local galaxies, and to point the 
way to similar analyses at higher redshifts. 
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M g log R ef f, g [kpc] 

Figure Al. R e g — M g (left) and n g — R c g (right) correlations. The small grey points are the datapoints in our sample, while the blue 
lines and shaded regions are the medians and 25 — 75 pe r cent scatter. The red line in the right plot is for the r-band data. Large grey 
points ar e for bright galaxies in [ Kormendv et al.l II2009T ). Black points with error bars are the _B- band results from T +09 (using the 
sample in IPrugniel fc Simienl <1996h l. while the cy an lines show th e resul ts using NYU-VAGC data iBlanton et al.ll2005ri , and the green 
points are the datapoints for the Virgo sample in iFerrarese et all lj2006T ). Suitable conversion factors are adopted to convert r-, B- or 
V-band data into g-band, and the data have been converted to our adopted cosmology. 



APPENDIX A: STRUCTURAL PARAMETERS 
AND MAGNITUDE 

Structural parameters of ETGs, such as the effective ra- 
dius, i? ft, are known to be significantly affected by the 
method (e.g. ID vs. 2D fitting of the light distribution) 
and profile shape (e. g. de Vaucouleurs vs. Sersic law) used 
to derive them (e.g. iKelson et al1 l2000r). While such differ- 
ences do not affect the FP relation of ETGs, because of the 
correlated variation of R a a and the mean surface bright- 
ness therein, (fj) e , they can affect other correlations among 
structural parameters and luminosity (mass) , and hence the 
correlations with DM fractions (densities). For instance, in 
Sec. I3.2l we have shown that the /dm correlations presented 
in the present study are somewhat steeper than those we 
previously deri ved in T+09, base d on a different sample 
of ETGs, from IPrugniel fc Simienl (| 19961 . hereafter PS96). 
We address this issue in Fig. IA1I where we compare our g- 
band luminosity-size relation to that of the PS96 B-band 
sample analyzed by T+09. For the sake of completeness, 
in the same plot we include our g-band Sersic index-size 
relati on, and compar e both relat ions to the g-ba n d dat a 
from IFerrarese etal] (120061 ) and iKormendv et al.l l|2009ri . 
Notice that the effective parameters of the PS96 sample are 
based on de Voucouleurs (rather than Sersic) galaxy pro- 
file fits, hence no comparison to our Sersic n-size relation 
is possible. We refer the reader to previous papers of the 
series (Papers I, II) for a detailed comparison of correla- 
tions among structural parameters from the SPIDER survey 
and SDSS pipeline. The present sample of ETGs exhibits 
a steeper luminosity-size relation than that of the T+09 
sample. This is explained by the fact that fitting a high-n 
galaxy with a de Vaucouleurs profile gives a systematically 
smaller R e g value, hence flattening the luminosity-size re- 
lation (e.g. iGraham fc Worlevll2008t l. Indeed, this explains 
why our /dm-RcS relation is steeper than that of T+09 
(see Fig. [5}, as a smaller R e ff encloses a smaller portion 
of the DM halo (implying lower /dm estimates). On the 



other hand, the correlations among structural parameters 
fo r the present sample a re in better agreement with those 
of IFerrarese et al.1 (|2006l ) and IKormendv et alj (|2009T l. We 
remark that since significant differences in the /dm correla- 
tions exist when analyzing different samples of ETGs (i.e. 
different sets of structural parameters), it is of paramount 
importance to compare different samples of galaxies (e.g. 
different environments) based on a homogeneous set of mea- 
surements, as the one used here. 



APPENDIX B: DYNAMICAL MASS 
PROCEDURE 

Our basic approach to estimating the DM content of 
ETGs is the same as in T+09 . For each galaxy, we 
first de-project the corresponding ISersid (|l968l ) light pro- 
file, assuming spherical symmetry. This gives the de- 
projected luminosity profile, j*(r), where r is the spa- 
tial distance to the centre of the galaxy. Using a Sersic 
model here has the main advantage of taking into ac- 
count the non-homolo gy of the stellar matter distribu- 
tion of ETGs (see, e.g., Caon. Capaccioli fc D'Onofridll993l 
and IPrugniel fc Simienll 19971 ) We also assume some simpli- 
fied functional form for the dynamical mass profile, M(r), 
i.e. either a const-M/L profile, M(r) — To L(r), or a SIS 
model, where M(r) oc ag is r (see Sec. [2}. Assuming spheri- 
cal symmetry and no rotation, we write the Jeans equation 
as: 



^ +2 m 2 oup 

dr r r 2 



(Bl) 



where [3 = 1 — a\ /o> is the anisotropy. Under the hypothesis 
of isotropy (i.e., /3 = 0), Eq. (|B 1 f) simplifies to 



a 2 r {r) = 



. GM j 



(B2) 
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In order to fit the given mass to the data, we project this 
equation in 2D, obtaining the line-of-sight velocity disper- 
sion: 

where 

p oo 

m = 2 L #F dr (B4) 

is the projected surface brightness profile. We integrate trios 
within a fixed aperture Ra p (i.e. the SDSS fibre aperture) 
to obtain the aperture velocity dispersion, <ta p : 

edp(flAp) = — ^- / ^ 2tv SI (S)aL(S)dS , (B5) 

where L(R) = J 2ttSI(S) dS is the luminosity within 
the projected radius R. To avoid lengthy calculations, we 
have adopted th e comp act fo rmulae for <tap c a lculat ed in 
Mamon fc Lokas] j2005af ) and iMamon fc Lokasl ^OOrj ) (see 



Mamon fc Loka s (2005b) for the anisotropic case). 

Finally, we fit the model <ta p to the observed a- spec , by 
varying the free parameters in Eq. (|B2|) (i.e. either crsis 
or To) until the desired matching is achieved. The result- 
ing best-fit mass profile provides the dynamical (spherical) 
mass-to-light ratio within R e a, Td yn (which coincides with 
To in the case of a const-M/L model), and mass M<j yn - Fi- 
nally, as we have shown in Sec. l3.ll we notice that a non null 
radial anisotropy (i.e. /3 > 0) decreases the Md yn , lowering 
the DM fractions (see text for further details). 



